Stochastic models for heterogeneous relaxation: Application to 
inhomogeneous optical lineshapes 

Gregor Diezemann, Gerald Hinze and Hans Sillescu 

Institut fur Physikalische Chemie, Universitat Mainz, Welderweg 15, 55099 Mainz, FRG 

§ ! Dynamic heterogeneity has often been modeled by assuming that a single-particle observable, 

fluctuating at a molecular scale, is influenced by its coupling to environmental variables fluc- 
tuating on a second, perhaps slower, time scale. Starting from the most simple Gaussian 
Markov process we model the exchange between 'slow' and 'fast' environments by treating the 
fluctuating single-particle variable as a projection from a higher- dimensional Markov pro- 
cess. The moments of the resulting stochastic process are calculated from the corresponding 
Master equations or Langevin equations, depending on the model. The calculations show 
the importance of the way to treat exchange processes. The resulting stochastic process is 
non-Markovian for all models. However, the deviations from a Gaussian behavior depend 
on the details of the models. A comparison of our results with other model treatments and 
^ ; experiments should provide further insight into the concept of dynamic heterogeneity. 
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I. Introduction 
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*^ ■ In spectroscopy inhomogeneous lineshapes are observed if individual molecules within a 
sample give rise to different spectral components. This can be observed even in isotropic 
systems if the relevant coupling depends on molecular variables, like orientations, that differ 
from each other. Another, more important source for inhomogeneous line broadening is 
provided by intermolecular couplings with the environment. These give rise to dynamic 
fluctuations and can thus be explored by lineshape studies. Classical examples are given 
by the motional or exchange narrowing of NMR spectra[|I|]. The Anderson- Weiss theory [^] 
of narrowing Gaussian lineshapes has been a model for many later treatments although 
it applies to homogeneous systems only. Furthermore, typical 'slow motion' lineshapes in 
NMR are far from being Gaussian. However, for our purposes, the Gaussian approximation 
provides a good starting point since it implies the well-known Ornstein-Uhlenbeck (OU) 
process^. This allows one to introduce non-Gaussian perturbations in a transparent way. 

A shortcoming of the OU process is that the two-time correlation function (2t-CF) of 
this Gaussian Markov process decays exponentially^. In complex systems like supercooled 
liquids, however, typically a non-exponential decay of 2t-CFs is observed. This forces us 
to consider non-Markovian models for the molecular observables. There are other fields in 
chemical physics where, starting from Gaussian fluctuations, non-Markovian effects have 
been incorporated in order to account for a more complex dynamics. Prominent examples 
are the problem of escape through a fluctuating bottleneck ||, ^ first order chemical 
reactions]?], || and single molecule spectroscopy 0. 

In our treatment, we closely follow the theory of solvation dynamics in time-resolved 
optical spectroscopy ^ , |TT| that has also previously been applied to treat heterogeneity in 
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terms of observed deviations from Markovian and Gaussian behavior |T2|, [13], [Uj. It should 
be noted that the average frequency, (u(t)), remains constant in essentially all dynamical 
NMR and ESR spectra but it changes due to time-dependent Stokes shifts in time-resolved 
optical spectra. This can result in time- dependent optical linewidths that have been related 
to dynamic heterogeneity [|12] , |i"5f . 

In the following Section, we briefly review the theory of time-dependent optical lineshapes 
in the form that will be used in the calculations. In Section III we treat the optical transi- 
tion frequency u(t) as a non-Markovian stochastic process, defined as the projection from a 
composite Markov process. For simplicity, we allow for two values of an additional environ- 
mental variable, e(t), to treat exchange between a 'slow' (e s ) and a 'fast' (e/) environment. 
Exchange in two (or multiple) state models has been treated a long time ago in relation 
with rotational [16|, [IT], IE, [L9| and translational[]20], |21| diffusion in complex systems. Here, 



we concentrate on the question regarding deviations from both, Gaussian and Markovian, 
behavior. The paper closes with a discussion of our results and conclusions in Section IV. 

II. Short survey of the theory of optical lineshapes 

This section consists of a brief review of the theory of time-dependent optical lineshapes 
in the inhomogeneous limit. In case of triplet state solvation all effects related to internal 
conversion, intersystem crossing, and the following electronic relaxation are assumed to be 
much faster than the dynamical features monitored experimentally. We mainly follow the 
general theoretical framework of refs.[|l(], [11], |14[| in the following. There, it has been shown 
that the time-dependent signal for a single chromophore immersed in a bath of solvent 
molecules can be written as: 

I(u,t)= du P(u,t\u )po(uj ) (1) 



Here, the transition frequencies u(t) depend on all relative distances and orientations be- 
tween the chromophore and the solvent molecules, which in turn depend on time. In Eq. (JJ) 
P(u,t\uo) denotes the conditional probability of finding a transition frequency uj at time 
t, given it was u at time t — 0. The initial probability po(u) is given by the emission 
spectrum at time t — 0, I(uj,0), because of P(u,0\uj ) = 5{oj — u ). Thus, observation 
of I(u,t), Eq. ([!]), as a function of time is interpreted as following the trace of the initial 
non-equilibrium emission spectrum towards the equilibrium steady-state emission spectrum 
I(u,oo) = p eg (uj) = P{uj, oo|u;o). Furthermore, it is assumed that the time-independent 
spectra are well approximated by Gaussians, 

where A = (u(0)), &q = ((u;(0) — A) 2 ) and cr 2 ^ = (u(oc) 2 ). The mean value (u(oo)) has 
been set to zero without loss of generality. Therefore, the overall red-shift is given by A. A 
Gaussian approximation for the conditional probability, 

P(u, t\oo ) = , exp -_L °WJ (3) 
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follows from a second order cumulant expansion with the Stokes shift correlation function 
given by C(t) = (u(t)w(Q)) / 'a 2 ^. Using Eq.(|3|) in Eq.(^J) for the time- dependent spectrum 
again gives a Gaussian for I(u, t) with a time- dependent mean and a second moment 

(u(t)) = AC(t) and (u, 2 (t)} = *l + (A 2 + a 2 - <£)C(t) 2 (4) 

yielding a time-dependent variance 

a\t) = (uj 2 (t)} - (co(t)} 2 = al + (a 2 - a 2 JC(t) 2 (5) 

This variance shows a monotonic behavior as a function of time in that it changes 
monotonously from cr 2 (0) = a 2 to a 2 (oo) = a 2 ^. 

It is precisely the deviations from the simple relation (|5]) observed in the triplet state 
solvation studies on supercooled liquids^, [TJ|] that has led to the interpretations of the re- 
sults in terms of heterogeneous dynamics. In supercooled liquids, the Stokes shift correlation 
function decays non-exponentially and may be written as 



C(t) = / dAp(A)e" 



(6) 



with relaxation rates A and a distribution function p(X). Performing a so-called inhomoge- 



neous cumulant expansion of the conditional probability, one finds WA 



P(u,t\uj ) = 
P w (uJ,t\u ) 



d\p(\)P(x)(u,t\u ) 
1 



: exp 



u e 



2a 2 (1 



(7) 



instead of Eq.(||). For the first moment, one again finds (w{t)) = AC(t), cf. Eq.(^), whereas 
the variance is given by 

° 2 (t) = al + (a 2 - al)C(t) 2 + A 2 [C(2t) - C(t) 2 } (8) 

This expression was found to give excellent agreement with the experimentally observed 
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It is seen that according to Eq.(|7]) the stochastic process u(t) can be viewed as originating 
from a superposition of Gaussian processes. This means that u(t) is modeled as a non- 
Markovian Gaussian stochastic process. In the next Section we discuss various exchange 
models, with particular emphasis on the deviations from Gaussian behavior. 



III. Calculation of frequency moments in exchange models 

The simplest stochastic model for a process like solvation is a Gaussian Markov process, 
i.e. the OU-process. If the transition frequency u(t) is modeled as an OU-process, one may 
consider the Langevin equation (LE): 

u(t) = - 7 u(t)+£(t) (9) 

where £(£) is a delta-correlated white noise, (£(£)) = 0, (£(£)£(£')) = 2c 2 7^(t — t') and 7 
denotes the damping rate. The LE, Eq.(^|), approximates the spectral diffusion process 
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as occuring on harmonic potential-energy surfaces. The transition frequency spreads in a 
restoring potential V = { r y/2)uj 2 . Alternatively, one can consider the corresponding Fokker- 
Planck equation (FPE) for the conditional probability: 



P(uj,t\uj ) = U FP (co)P(uj,t\uj ) with IT 



FP, 



10) 
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(10) 



The resulting expression for P(u>, t|co>o) is given by Eq.(|3]) with C(t) =exp (— jt). Therefore, 
in order to have a model which yields a Stokes shift correlation function decaying non- 
exponentially according to Eq.(^), one has to consider u(t) as a non-Markovian stochastic 
process. In the present paper, this is achieved by considering exchange models. As already 
indicated in the introduction, we restrict ourselves to two-state models, because these allow 
for analytic solutions. The idea underlying such models can most easily be summarized in the 
following way. We consider the solute to experience various surroundings characterized by an 
environmental variable e^, i G {s, /}. Thus, in a given environment, the spectral diffusion of 
the optical transition frequencies u(t) occurs with a decay constant 7« = 7(ej). However, the 
environments are not viewed to be static, but rather fluctuate in the course of time. These 
environmental fluctuations are characterized by exchange rates K^f. = K,(ei\ek). In order to 
further simplify the calculation, we take the exchange rates to be equal, k s j — Kf jS — j x , 
implying Pi q = p eq {u) = 1/2. The generalization to a larger number of states is trivial. Also 
other models, like energy landscape models|23|], can be considered |T4 



Even though we restrict ourselves to simple two-state models, we have to face the am- 
biguity concerning the way of how to incorporate the exchange between the fast and slow 
environments into the model. In the following, we consider two different, but in some way 
strongly related ways of treating the rate exchange. 



A. Master equation models 

One class of models allowing to treat rate exchange processes can be defined by supple- 
menting the FPEs in the slow and fast environments with some terms representing the 
exchange process. This way one obtains rate equations for the conditional probability 
P ijk (u,t\uo) = P(u,€i,t\u} ,€k) for finding the pair {uj,€i} at time t, given it had values 
{u , tfc} at t = 0. Still, we have to make some additional asumptions about what happens 
with the value of to during the exchange process. Two limiting choices for this are the 
following [[17], [|{|. Either one assumes that u does not change at all during a ^ — > e k transi- 
tion. Alternatively, one allows u to randomize completely during such a transition, meaning 
that after the transition u can take any value according to its statistical weight p eq (uj). Let 
us start with the first choice. This is the same assumption that has often been used in 



context of models for rotational and translational motions in supercooled liquids [HO. 21] and 



also in treatments of spectral diffusion] 13 



We abbreviate the FP-operators in the fast and in the slow environments by Tlf p (u). 
These operators are defined in Eq.flTOD, but with ^ replacing 7. We have to solve the coupled 
FPEs: 

Pij(u,t\wo) = Tl[ p (uj)P i! i(uj,t\uJo) - 7^ [P»,/(a;,*|a;o) - P ki i(u,t\u )] ; i^k (11) 

Due to the coupling of the FPEs, the resulting Pi t i(u,t\uo) are not Gaussians. Furthermore, 
an analytical solution of Eq.(|ll|) is not feasible. In order to calculate the spectral moments 
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or any other observable, an expansion in terms of the eigenfunctions of the FP operator|M| is 
utilized, P^{u, t\uj ) = e~ w2 J2^=o N n H n (u)H n (u )G < f' f !(t), where H n (z) are Hermite polyno- 
mials, u 2 = u 2 /2a 2 <1 and N' 1 = 2 n n\(2na 2 <> ) 1 / 2 . Using this expression in the FPEs, Eq. QXTj) , 
results in the following equations for the Green's functions: 



d_ 
dt 



G&\t) = - (»7i) G%\t) - 7, GS'(t) - where k?i (12 



which can be solved easily. In the following, we assume o"o = (J^ for simplicity. For the 
lowest-order spectral moments one finds, after performing the cu-integrations: 

(u(t)) = AC 1 (t) 

(u, 2 (t)) = al + A*C 2 (t) (13) 

Furthermore, one has (u 3 (t)} = 3A a^dit) + A 3 C 3 (t) and (u\t)) = 3^ + 6a^A 2 C 2 (t) + 
A 4 C*4(t). In principle, all moments are accessible this way, although the calculations be- 
come rather tedious. It is important to notice at this point that the n th moment (u n (t)) is 
determined by the functions 

Cn(t) = ~ E G$(t) = pW exp [- 7 W*] + P f exp [-yft] (14) 

i,k 

Here, the pf 1 ^ and 7}^ are determined by the solution of Eq. (|i~2|) : 

P$ = i (l ± ) 7$ = ~ (27, + n( 7/ + Is) T Z n ) (15) 
with Z n = V / (n 2 ( 7/ - 7s ) 2 + 4 7 2) 

In Fig. la we plotted the Stokes shift correlation function Ci(t) = C(i) versus time. For 
the parameters chosen it is clearly seen that mainly the long time decay is affected by a 
finite exchange rate j x . Fig. lb shows the time-dependent part of the variance, a 2 c (t) = 
[u 2 {t) — a^j/A 2 versus time for the same parameters. Additionally shown is the result of 
the Gaussian approximation, C\(2t) — Ci(t) 2 , cf. Eq.(S). Remember, that this expression 



is in excellent agreement with the available experimental data[H|. It is clearly seen that 
a 2 sc (t) exceeds this expression, in particular in the time range, where both functions exhibit 
their maximum. It is exactly this excess observed for this particular model that previously 
led to the incorrect conclusion that the experimental data are only compatible with a long- 
lived dynamic heterogeneity jl3| . In the present model the differences between a 2 c (t) and 
Ci(2t) — Ci(t) 2 originate solely from the fact, that the stochastic process u(t) is not only 
non-Markovian, but also non-Gaussian. 

From the expressions given in Eq. flT5|) for n ^> 1 one approximately has 7^™j ~ 7^ + nj s / f 

and p™\ ~ 1/2 ± j x / (n(7/ — 7s )). Therefore, using the parameters of Fig.l, one finds that 

even for n = 100 the exchange rate 7^ still contributes on the order of 30% to 7 ^ 100 ' 1 . In Fig. 2 
we show the quantities ( n (t) = C n {t) — C\{nt). These give a measure for the deviations from 
a Gaussian behavior because in that case one has ( n {t) = for all n > 0. This plot nicely 
illustrates the fact just discussed. With increasing n the maximum ( n (t) is shifted towards 
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shorter times t max ^ n , Fig. 2a. However, one can see that the maximum of Cn(t) roughly occurs 
at nxt max ^ n ~7^T 1 , as shown in Fig.2b. This shows that not only the first few moments are 
affected by the deviations from Gaussian behavior in this model, but that these deviations 
even increase with n. 

In the above treatment we considered the case in which the transition frequencies u(t) do 
not change during a £j — > transition. We already discussed that another extreme scenario 
would be defined by allowing u(t) to randomize during such a transition. We just mention 
that for this model the eigenvalues for n = remain the same as in the case considered above. 
For finite n, however, one has p^X = 1/2 and 7^ = jx+n^fs/f- Thus, the qualitative features 
are not very different in this case. In particular, in these exchange models the stochastic 
process uj (t) is non-Markovian and non-Gaussian. 

B. Langevin equation models 

Instead of starting from the coupled FPEs, Eq. (|TT|) , one can also define an exchange model 
by considering coupled Langevin equations. In such models one associates a value 0Ji{t) with 
either of the states i = s, f and introduces an additional exchange term, cf. Eq.(|9]): 

Ui(t) = - (ji + J x ) LOi{t) + J x UJ k (t) + £i(t) (16) 

where again i 7^ k and is a delta-correlated white noise, meaning that (£«(£)) = 
and (£,i{t)!;k{t')) = r ii / c 5(t — t'). Note, that here we do not have different choices for the 
incorporation of exchange effects, since by definition the LE model involves diffusive motions 
only. Therefore, in such a model u changes by an infinitesimal amount during an 6j — > 
transition. The solution of Eq.flT6"D results in two-dimensional Gaussians for the conditional 
probabilities. In order to show how the results are related to the ones obtained for the FPEs, 
we proceed in the following way. We introduce the eigenvalues X a and the corresponding 
eigenvectors \a) of the matrix W defined in Eq. ([T6|) . In such a notation we have W^k = 
(i|W|fc), W s , s = -7, - 7-,., W f j = -7/ - 7s, W sJ = W ftS = j x and u^t) = (i\uj(t)). 
Correspondingly, using J2i = J2 a \ a )i a \ = 1 an d the orthonormality of the vectors, we 
have cjj(i) = J^aiA^^ait) and Wi^ = Ha{A°) . This yields the decoupled LEs: 

u a (t) = \ a u Q (t) + £ a (t) (17) 

from which it is imediately evident that the solutions for the corresponding P a (u>, i|o?o)0 
are Gaussians determined by the correlation functions C a {t) = (L} a (t)u a (0)). In Eq . (|P7D , the 
projections of the noise give {£ a (t)£ a (t')) = T a 5(t — t'). Choosing T^k = — 2<T 00 W / j i fc yields 
r a = — 2X a . For the moments one finds: 

Mt)) = MO))^ and K(t)) = (^(0))e 2 ^ + ^(l-e 2 ^) (18) 

and the correlation function is given by C a (t) = (^(0)) e Xat . For an arbitrary observable the 
expectation value is given by (A(t)) = (1/2) J2a(A a (t)) , and thus 

(u{t)) = A d(t) and (cu 2 (t)) = + A 2 d(2t) (19) 

where the C n {t) are defined in Eq . fll4]) . Therefore, the time-dependent variance in this model 
is given by a 2 (t) = a 2 QO + A 2 [d(2t) - Ci(t) 2 ], cf. Eq.©. In Fig.lb, C x {2t) - d{t) 2 is shown 
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as the dashed line. The mathematical reason for the behavior of the moments, Eq.flI5|), lies 
in the fact that the matrix W defined in Eq.(^) is identical to the one for Gfl(t) defined 
in Eq.([T^). This means that the eigenvalues A a are just given by X s /f = —%h according 
to Eq. ([T5|) . The conclusion from this observation is that the Langevin model considered 
here is nothing but the Gaussian version of the master equation model studied above. A 
simple calculation shows that the LEs, Eq . flTJD , considered here corresponds to the coupled 
FPEs P iyl (uj 1 t\ujQ) = Y, k tlf k p {uj)P Kl (uj,t\uj ). Here, fl[ k p {uj) is obtained from Eq.fllOD by the 
replacement 7— W^fc). The important point is that in this model uj{t) is a non-Markovian 
Gaussian process. 



IV. Discussion and Conclusions 

In the preceeding section we have analyzed simple exchange models that allow one to account 
for fluctuations in dynamic heterogeneous systems. In the particular example of solvation 
dynamics we have focussed on the calculation of the lowest-order moments of the transition 
frequency distributions. Despite its simplicity this example allows to discuss various aspects 
of dynamic heterogeneity. In particular, without exchange the transition frequency u(t) con- 
stitutes a Gaussian Markov process, the OU process. This is among the simplest stochastic 
processes one can study 0. 

We have investigated different variants of the models to incorporate the way in which 
exchange processes change the temporal evolution of the conditional probability. These 
exchange processes are modeled as transitions between different environments, characterized 
by different values of an additional stochastic process e(t). In the master equation models 
one has to make specific assumptions about the effect of exchange on the stochastic process 
u>(t). We discussed only two extreme choices. One choice, which also has been used in 
many earlier investigations, corresponds to the assumption that the processes u{t) and e(t) 
are uncorrelated. This means that u(t) is not at all affected by a change in e(t), i.e. by 
exchange. We have shown that the equations determining the moments (u n (t)) can then be 
solved analytically. This allows to discuss their behavior in some detail. A very important 
finding is that in this case the process u(t) is neither Gaussian nor Markovian for finite 
exchange rates j x . In contrast, when studying a Langevin model, we found that in that 
case uj(t) also is non-Markovian but Gaussian. Therefore, using this model yields excellent 



agreement with the experimental data|15[ for arbitrary choices of the exchange rates. This, 
however, is not astonishing because according to Eq.(0) any model that yields a description 
of uj(t) in terms of a non-Markovian Gaussian stochastic process will give results compatible 
with the experimental data. 

The most important point in the discussion of the lifetime of the dynamic heterogeneities 
appears to be the question as to how to quantify this lifetime. It has been stressed earlier 



that this cannot be achieved by any two-time correlation functionp5||. This fact stimulated 
the development of NMR techniques to study higher-order time correlation functions |2(J. A 
comparison of the conditional probabilities P(u, t\uo) with two-dimensional NMR spectra |27| 
shows that these are identically the same if u is reinterpreted as the spin-precession fre- 
quency. The mentioned time dependence of (u)(t)) that is absent in NMR does not enter 
the discussion here, because this time dependence does not alter P(ou,t\uj ). Thus, even if 
the full conditional probability would be known, one could not obtain information about 
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the heterogeneity lifetime. However, knowledge of higher moments would allow to better 
specify the detailed 'geometric' properties of the spectral diffusion process. In particular, 
one could get information about possible deviations from diffusive behavior, which was our 
starting point in this paper. It is exactly this issue which can be treated very efficiently by 
2d-NMR methods. One can clearly distinguish between rotational diffusion and rotational 



jump models p7 



The comparison between the spectral diffusion process and the reorient at ional motion of 
molecules in supercooled liquids further shows the different impact of exchange on various 
models. Rotational diffusion (or also rotations by finite jumps) represents a Markov process. 
This process, however, is not Gaussian. Therefore, in this example it hardly matters in which 
way exchange processes are incorporated. By contrast, in the example of spectral diffusion, 
different ways of treating exchange give different results concerning the properties of the 
resulting stochastic process. As discussed above, using master equation models usually will 
destroy the Gaussian properties of oj(t). Therefore, in this case one has a larger variety of 
modifications that can be achieved by considering various exchange models. 

In conclusion, we have investigated different scenarios for exchange processes in dynam- 
ically heterogeneous systems. Starting from a Gaussian Markov process, we have shown 
that different assumptions concerning the coupling of u>(t) to the exchange processes yield 
different results. In particular, if the values of u change continuously during an exchange 
process the resulting stochastic process is non-Markovian but still Gaussian. Such a 'diffu- 
sive' change of u(t) during exchange of course means that the effect of 'exchange' (7^) is 
very similar to the one of 'relaxation' (7^//). In contrast, if exchange does not at all affect 
the values of u or if it causes them to change discontinuously also the Gaussian properties 
of the resulting stochastic process is lost. In any case, deviations from a Gaussian behavior 
are not good candidates for the determination of heterogeneity lifetimes. Still, they may be 
helpful in discriminating among various variants of exchange models. 



Acknowledgement 

We are grateful to R. Bohmer for fruitful discussions and to the Deutsche Forschungsge- 
meinschaft for continuing support of our research on dynamic heterogeneities via the Son- 
derforschungsbereich 262. 



8 



References 

[1] A. Abragam; The Principles of Nuclar Magnetism, Clarendon, Oxford (1961) 

[2] P.W. Anderson and P.R. Weiss; Rev. Mod. Phys. 25 269 (1953) 

[3] N.G. van Kampen: Stochastic Processes in Physics and Chemistry, North- Holland, 
Amsterdam, New York, Oxford (1981) 

[4] R. Zwanzig; J. Chem. Phys. 97 3587 (1992) 

[5] J. Wang and PG. Wolynes; Chem. Phys. Lett. 212 427 (1993) 

[6] D.J. Bicout and A. Szabo; J. Chem. Phys. 108 5491 (1998) 

[7] N. Agmon and J.J. Hopfield; J. Chem. Phys. 78 6947 (1983) 

[8] J. Wang and PC Wolynes; Chem. Phys. 180 141 (1994) 

[9] J. Wang and PC Wolynes; Phys. Rev. Lett. 74 4317 (1995) 
[10] R.C Loring, Y.J. Yan and S. Mukamel; J. Chem. Phys. 87 5840 (1987) 
[11] M.D. Stephens, J.C Saven und J.L. Skinner; J. Chem. Phys. 106 2129 (1997) 
[12] R. Richert; J. Chem. Phys. 114 7471 (2001) 
[13] R. Richert; J. Chem. Phys. 115 1429 (2001) 
[14] G. Diezemann; manuscript submitted to J. Chem. Phys. 
[15] HWendt and R.Richert; Phys. Rev. E 61 1722 (2000) 
[16] J.E. Anderson and R. Ullman; J. Chem. Phys. 47 2178 (1967) 
[17] D. Becker and H. Pfeifer; Ann. Phys. 16 22 (1965) 
[18] H. Sillescu; J. Chem. Phys. 54 2110 (1971) 
[19] H. Sillescu; J. Chem. Phys. 104 4877 (1996) 
[20] J. Karger; Ann. Physik, 7. Folge 24 1 (1969) 

[21] I. Chang and H. Sillescu; J. Phys. Chem. 101 8794 (1997); E. Bartsch, T. Jahr, A. 
Veniaminov and H. Sillescu; J. Phys. IV, France 10 289 (2000) 

[22] R. Richert; J. Chem. Phys. 113 8404 (2000) 

[23] G. Diezemann; J. Chem. Phys. 107 10112 (1997) 

[24] H. Risken: The Fokker-Planck Equation, Springer, Berlin (1984) 



9 



[25] R. Bohmer, R.V. Chamberlin, G. Diezemann, B.Geil, A. Heuer, G. Hinze, S.C. Kue- 
bler, R. Richert, B. Schiener, H. Sillescu, H.W. Spiess and M. Wilhelm; J. Non-Cryst. 
Solids 235-237 1 (1998) 

[26] K. Schmidt-Rohr and H.W. Spiess; Phys. Rev. Lett. 66 3020 (1991) 

[27] K. Schmidt-Rohr and H.W. Spiess; Multidimensional Solid-State NMR, Academic 
Press, London 1994 

Figure captions 

Fig.l : The first and second moments versus time for the master equation model. 

a: The Stokes shift correlation function Ci(t) according to Eq.flTJP versus time for 
7s = 0.1s _1 , 7/ = 10.0s _1 , 7 x = 0.33s _1 (full line) and 7^ = (dotted line), 
b: The time-dependent part of variance, c 2 sc {t) = [o~ 2 {t) — er^J/A 2 , versus time (full 
line) for 7 X = 0.33s -1 . Also shown is C\(2t) — Ci(t) 2 as dashed line. This corresponds 
to the time- dependent part of the variance in the Langevin equation model. The 7 S 
and 7/ are the same as in (a). 

Fig. 2 : The quantities ( n {t) = C n (t) — C\(nt) versus time for the same parameters as in 
Fig.l ( r y x = 0.33s~ 1 ). a: ^ 2 (^) : full line, Cio(^) : dashed line, C2o(^) : dot-dashed line and 
Cioo(^) : dotted line. 

b: The same ( n (t) as in (a) are plotted versus the scaled time (nt). The arrow indicates 
7^ = 3.05. 
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